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We have demonstrated transport of 9 Be + ions through a 2D Paul-trap array that incorporates 
an X-junction, while maintaining the ions near the motional ground-state of the confining potential 
well. We expand on the first report of the experiment [l] , including a detailed discussion of how the 
transport potentials were calculated. Two main mechanisms that caused motional excitation during 
transport are explained, along with the methods used to mitigate such excitation. We reduced the 
motional excitation below the results in Ref. [T] by a factor of approximately 50. The effect of a 
mu-metal shield on qubit coherence is also reported. Finally, we examined a method for exchanging 
energy between multiple motional modes on the few-quanta level, which could be useful for cooling 
motional modes without directly accessing the modes with lasers. These results establish how 
trapped ions can be transported in a large-scale quantum processor with high fidelity. 
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I. INTRODUCTION 

The reliable transport of quantum information will en- 
able operations between any arbitrarily selected qubits 
in a quantum processor and is essential to realize effi- 
cient, large-scale quantum information processing (QIP). 
Trapped ions are a promising system in which to study 
QIP [2H3], and several approaches to achieving reliable 
information transport have been proposed 042] ■ In most 
demonstrated entangling gate operations that use ions, 
qubits stored in the internal atomic states of ions are 
entangled by coupling the internal states with a single 
shared motional mode through a laser-induced interac- 
tion pHU HQ]- However, as the number of ions grows 
large (> 10), it becomes difficult to isolate a single mode 
during gate operations [51 lllj . One way around this issue 
is to distribute the ions over an array of harmonic poten- 
tials, where the number of ions in each trapping potential 
can remain small. The potentials can be adjusted tem- 
porally to transport the ions throughout the array and 
combine selected ions into a particular harmonic poten- 
tial. Once combined, gate operations can be performed 
on the selected ions by use of a local shared mode of 
motion (BJE]. 

Initial demonstrations of such distributed architectures 
have incorporated simple linear arrays [12}jl6j, where all 
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ions are confined in potential minima on a line along 
an axis of the trap. The order of ions within the lin- 
ear array can even be changed [T7]. However, multidi- 
mensional arrays [6j |7] provide the greatest flexibility in 
ion-trap processor architectures, and permit more effi- 
cient reordering of ion strings for deterministic gate op- 
erations. The key technical element that must be realized 
towards this end is the two-dimensional junction, which 
consists of multiple intersecting linear arrays. The po- 
tentials in a junction are more complicated than those 
in a linear array, making transport through a junction 
challenging. Since the fidelity of the gates is highest if 
the ions are near their motional ground state, it is im- 
portant that transport through such arrays be performed 
reliably and with minimal excitation of the ion's motion 
in its local trapping potential. If multiple transports are 
needed, each transport should contribute well under a 
single quantum of motional excitation, though sympa- 
thetic cooling can be used to remove excess motional en- 
ergy, at the cost of increased experiment duration (and 
accompanying decoherence) [7]. For simple linear arrays, 
reliable transport with little motional excitation has been 
demonstrated pi fT3l ITS] . 

To date, transport through a T-junction [TH], an X- 
j unction [1] and surface-electrode Y-j unctions [El [20] 
have been demonstrated. However, such transport has 
not yet been demonstrated with sufficiently-low motional 
excitation (at or below a single quantum). Using the 
apparatus in Ref. [I], we have now realized highly reli- 
able transport through an X-junction with excitation of 
less than one quantum of motion per transport, a de- 
crease of approximately 50 compared to the results in 
Ref. pQ. This has allowed us to observe a process where 
energy can be exchanged between motional modes in cer- 
tain situations, and demonstrates motional control over 
the ions at the single-quantum level. The paper is orga- 
nized as follows: we begin in Sec. [H] with a description 
of the X-junction trap array used for transport. Sec- 
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tion |III| lays out the procedure for calculating the time- 
dependent trapping potentials that transport the ion. A 
description of the basic transport experiment is given in 
Sec. IIVI Section |V] covers the various mechanisms that 
excite the ion's motion during transport, as well as the fil- 
tering techniques used to mitigate those excitations. This 
understanding of the noise sources, and the subsequent 
improved filtering techniques, allowed the reduction in 
motional excitation relative to Ref. pQ. To mitigate the 
effects of magnetic field fluctuations on qubit decoher- 
ence, a mu-metal shield and field-coil current stabiliza- 
tion were used, which is explained in Sec. |VI| Finally, in 



Sec. VII we discuss a procedure for swapping motional 
energy between motional modes at the center of the junc- 
tion array. This swapping process can potentially be used 
to laser-cool multiple modes of motion without the need 
for a direct interaction between the cooling laser and ev- 
ery motional mode. 



II. X-JUNCTION ARRAY 
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FIG. 1: A cross sectional view (not to scale) of the five-wafer 
stack, in the x, y plane at the experiment zone (£). Each wafer 
had a channel cut through it to define the electrode structure 
and to provide a path for laser beams to pass through the 
wafer stack. The top and bottom wafers provided the confin- 
ing potential; the ions were trapped between these electrodes 
as indicated. The RC low-pass filters were surface-mounted 
to the filter board with gold ribbon attached by use of resis- 
tive welds. The bias wafer was a single electrode used to null 
stray electric fields along y. Gold (represented in yellow) was 
coated on the top side of the trap wafers and wrapped around 
to the bottom side, and vice-versa for the bias wafer. Gold 
wire bonds connected traces on the trap wafers to traces on 
the filter board. 



The X-j unction array was based on the design of pre- 
vious two-layer linear RF Paul traps [H EH HQ- The 
current trap consisted of a stack of five high-purity alu- 
mina (99.6% AI2O3) wafers clamped together (Fig. I]) 
with screws (visible in Fig. [5]). The trap electrodes 
resided in the 'top' and 'bottom' wafers. These wafers 
(125 /im thick) were laser machined to cut out 'main 
channels' through the wafers, with opposite sides of the 
channel forming rf and control electrodes. Slits, nomi- 
nally perpendicular to the main channel axes, separated 
the control electrode side of the channel into a series of 
cantilevered structures to produce separate control elec- 
trodes. Electrodes were formed onto the AI2O3 by evapo- 
rating through a shadow-mask a 30 nm titanium adhesion 
layer followed by 0.5 /jm of gold, then overcoating with 
3 /xm of electroplated gold. Care was taken to coat all 
sides of each cantilevered structure to minimize exposed 
dielectric that could otherwise charge and shift the po- 
tential minima in an uncontrolled way. 

A spacer wafer provided a separation of 250 /im be- 
tween the two trap electrode layers. These three wafers 
sat atop a 500 /im thick 'filter board', upon which in- 
vacuum RC filtering components were mounted. The 
'bias wafer' resembled the 'top' and 'bottom' wafers 
but with a single continuous control electrode extend- 
ing along all sides of the main channels. The bias wafer 
sat below the filter board and was used to compensate 
stray electric fields along y. 

Gauge pins were used to help align the wafers dur- 
ing assembly. A misalignment error of approximately 
0.22° was measured between the z axes of the two elec- 
trode wafers, and this error was included in the com- 
puter model of the trap used to determine the appropri- 
ate transport potentials. 

The electrode layout of the array is depicted in Fig. [3] 
and consisted of 46 control electrodes that produced 
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FIG. 2: Top view of the filter board and trap wafers. The 
filter board fills the entire image, while the top wafer is the 
rotated square visible on the right of the image. Cap screws, 
visible in two corners of the top wafer, held the wafer stack 
together. Wire bonds connected the filter board traces to 
the top and bottom trap wafers. Surface-mount resistive and 
capacitive elements on the filter board provided filtering for 
the control potentials (see Fig. 10 1. 



18 possible trapping zones. The experiment zone, £ , was 
chosen as the zone where the ions interacted with lasers 
for cooling and qubit operations. In addition to £ , zones 
J 7 , V, and C (at the center of the junction) composed 
the four destinations of the transport protocols. The fi- 
nal zone of interest was the load zone, C, where the ions 
were initially trapped. 

The trap dimensions were similar to those in Refs. [13l 
|2T] . The width of the channel between the rf and control 
electrodes was 200 /im, except near C, where it increased 
to 300 /im to increase the volume of the loading zone and, 



3 




FIG. 3: (a) Cross-sectional view of the two layers of electrodes 
in the X-junction array, (b) Top view of the electrode layout, 
with the rf electrodes indicated, and all other (control) elec- 
trodes held at rf ground. The bottom trap wafer, which sat 
below these electrodes, had a nearly identical set of electrodes 
but with rf and control electrodes exchanged across the main 
channel. Ions were trapped in the main channels between the 
rf and control electrodes. Forty-six control electrodes (some 
of which are numbered for reference) supported 18 different 
trapping zones. The load zone (£), the main experiment zone 
(£), the vertical zone (V), the horizontal zone (J-) and the 
center of the junction (C) are labeled, (c) Schematic of the rf 
bridges from an oblique angle (not to scale). 



with it, the loading probability. Most control electrodes 
extended 200 /im along the trap axis, but those nearest 
to the junction were 100 /im wide to ensure sufficient 
control when ions were transported in this region. 

At C, two main channels crossed to form an X-junction, 
and two rf bridges connected the rf electrodes on opposite 
sides of that junction (one on the top trap wafer and one 
on the bottom). Without such bridges, the array would 
not have provided harmonic three-dimensional confine- 
ment at the center of the junction [2"2l [2"3"] . The widths 
of the bridges were 70 /im, though the trapping potential 
was not strongly dependent on this dimension. 

These bridges introduced four axial pseudopotential 
barriers, one in each of the entrances to the junction 
(along ±x and ±£). Figure [4] shows the two simu- 
lated pseudopotential barriers along the z legs in the 
X-junction array going toward £ and T (the asymme- 
try was due to the trap misalignment mentioned above). 
The height of these barriers was a significant fraction 
of the transverse pseudopotential trapping depth and 
was approximately 0.3 eV for 9 Bc + , with rf poten- 
tial of V T f fa 200 V (peak amplitude) and frequency 
r2 r f s» 27T x 83 MHz. At the apex of the barriers, just 
outside the center of the junction, the pseudopotential 
was anti-confining in the axial direction but still har- 
monically confining in the two radial directions. It was 
possible to use the control electrodes to overwhelm this 
anti-confinement and produce a 3D harmonic confining 
potential at all points along the axis of the array. 

Zone £ was positioned far (880 /im) from the junction 
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FIG. 4: Simulated pseudopotential barriers along the z axis 
produced by the rf bridges in the X-junction, with being 
the junction's center. Here, we assumed V r i ~ 200 V and 
fi r f ~ 2n x 83 MHz. The asymmetry between the two barriers 
was due to a slight misalignment of the trap wafers. 




FIG. 5: A Be oven was positioned out of the plane of this fig- 
ure in the positive y direction (above the trap) and could be 
heated to produce a flux of neutral Be. This Be would then 
travel down onto the trap, with a portion of the flux pass- 
ing into trap's main channel at the load zone. Copropagating 
photoionization and Doppler-cooling laser beams intersected 
the Be in the load zone at 45° to the xz plane of the page and 
parallel to — y + z. An 'L'-shaped oven barrier obscured the 
line-of-sight between the oven flux and the zones used during 
the transport experiments to prevent neutral Be from accu- 
mulating on the surfaces of the electrodes near the junction. 
This barrier was positioned just above the trap electrodes, 
extending 1.6 cm along y out of the plane of the page. Addi- 
tional laser access was available for beams passing through £ 
(at 45° to the xz plane) allowing for cooling, detection, and 
gate operations at £ . 

to reduce the residual slope of the pseudopotential barrier 
in this zone. The amplitude of the axial pseudopotential 
at £ was estimated, by use of computer models, to be 
2.9 x 10~ 5 eV with a 8.7 x 10~ 8 eV//xm axial gradient, 
which would give rise to an axial 'micromotion' amplitude 
of 47 nm at the drive amplitude specified above. 

Ions were loaded into the array from a flux of neutral 
Be that passed through L and was photoionized with a 
mode-locked laser that after two stages of doubling pro- 
duced 235 nm resonant with the S-to-P transition of neu- 
tral Be. To help prevent buildup of neutral Be from the 
beam in other regions of the array, C was located suffi- 
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ciently far from £. In addition C was displaced along x 
from £, by use of two 135° bends in the main channel, 
to allow an 'L'-shaped stainless-steel shield to be placed 
0.5 mm above the trap wafers, preventing neutral Be from 
striking the experiment zone while allowing laser access, 
as shown in Fig. [5] Transporting ions through such 135° 
bends is relatively straightforward, and we were able to 
easily transport ions between C and £ . 

Whenever an ion was lost, a new ion was loaded into 
C and immediately transferred to £ . It was also possible 
to use zone £ as a reservoir zone, where extra ions were 
loaded and held in reserve until needed to replace ions 
lost in the experiment region. This allowed the loading 
process to be performed less often, which avoided heating 
the neutral Be oven and the concomitant degradation of 
the vacuum. Potentially many ions could be simultane- 
ously stored in such a zone, though we only stored a small 
number and did not regularly make use of this feature of 
the trap. By enabling a better vacuum, a reservoir can 
significantly increase ion lifetime. In this scenario, it is 
desirable to maintain Doppler cooling in C to extend the 
ion lifetime. 



III. TRANSPORT POTENTIALS 

The first demonstrations of ion transport in a multi- 
zone trap involved moving an ion along a linear array |12j . 
A protocol where two ions were placed in a single trap- 
ping well and separated into two wells or combined from 
two wells into a single well was also demonstrated p~2|. 113] . 
Since then, transport through linear arrays has been ex- 
tended to other contexts [21 [THl HZ] i including transport 
through a junctions [U fTM2Tj] and switching of ion or- 
der dang. 

Here we outline the process used to calculate the time- 
series of control potentials, or 'waveforms', used to trans- 
port ions through the X-junction. This same basic proce- 
dure would be generally applicable to many ion-transport 
situations. The goal was to move ions quickly, over long 
distances, while maintaining low excitation of the ion's 
secular motion in its local potential, and traversing non- 
trivial potential landscapes such as those near junctions. 
Ideally, the ion should move along the axial direction of 
the array while remaining at the transverse pseudopoten- 
tial minimum. The control electrodes were used to create 
an overall harmonic trapping well whose minimum moved 
along this desired trajectory. The procedure for deter- 
mining waveforms can be broken down into four steps: 
modeling the trap, determining the constraints, solving 
for the appropriate potentials, and assigning the time de- 
pendence of potentials. 

An electrostatic model of the trap was constructed by 
use of boundary element method (BEM) software [Ml 
[25] . For each of ./V electrodes, the model was run once, 
applying 1 V to the n th given electrode while grounding 
all other electrodes. The potential resulting from each of 
these voltage configurations, (j> n (r), was extracted (in the 



form of a 5 /xm grid) in the region through which the ion 
would pass. These individual electrode potentials could 
then be weighted by the actual voltage applied to the 
electrode, V n , and summed to find the total potential: 



iV 



a) 



71=1 



Here we have included the contribution of the rf pseu- 
dopotential ps . The pseudopotential was found by first 
modeling the rf potential </> r f as if it were a static poten- 
tial at 1 V. Then an additional step was used to convert 
the rf potential into a time-independent pseudopotential 
by use of 



^p S (r) = 



g 



V rf V<Mr) 



(2) 



where V T { was the peak voltage applied to the rf electrode, 
and q and m are the charge and mass of 9 Be + , respec- 
tively. (Throughout this section, all <f> potentials (includ- 
ing 4> vs ) are reported as electric potentials (in units of V) 
and not energy potentials (units of eV); these are related 
by a factor of q.) 

The waveform was built up from a string of individ- 
ual solutions, where each solution satisfied a set of con- 
straints on the trapping potential centered at a certain 
position. These constraints are defined below, but relate 
to defining the secular frequencies and orientation of the 
principle axes of the potential. By advancing that posi- 
tion by 5 /im along the intended ion trajectory for each 
subsequent solution, the series of potential steps was cre- 
ated that moved the potential well along the sequence 
of positions. In theory, the constraints can be set to 
completely define a harmonic potential localized at the 
desired position, while also constraining the three secu- 
lar frequencies and the orientation of the principal axes. 
This would imply nine constraints, which we assume for 
now, though below we will relax some of these constraints 
when solving for the experiment waveforms. 

To produce a trapping potential, $(r), with a mini- 
mum at ro = (xo,yo, zq), we enforce 



V$(r ) =0, 

where = is used to mean 'constrained to be true'. 
The Hessian matrix, 



(3) 
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a 2 
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*(ro), 



(4) 



can be used to extract the remaining six parameters 
of the harmonic potential: the eigenvalues A^ of H(r ) 
are related to the secular frequencies, A^ — mwf and 
the eigenvectors point along the principal axes. By 
completely constraining the Hessian, we constrain these 
quantities. Note that the Hessian is symmetric (H = 
H 1 "), and has only six independent entries. 
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It is most convenient to evaluate the Hessian in the 
frame of the desired principal axes, (x',y',z'), in which 
case the Hessian constraint equation simplifies to 



and 



H(t ) = m 





,2 

o 



(5) 



where diagonal entries constrain the desired secular fre- 
quencies (u) x > , ujyt , oj z i ) , and the off-diagonal entries con- 
strain the principal axes to point along (V, y' , z'). If the 
Hessian is evaluated in a different basis, the right-hand 
side of Eq. ([5| will not be diagonal, and the frequency and 
axis constraints are mixed. Nonetheless, an appropriate 
choice for the right-hand side can still be made in that 
case [26]. From here on, we assume [x, y, z) = (x 1 , y', z'). 

In the interest of compact nomenclature, it is conve- 
nient to define several column vectors: 



V = [ 1 Vi v 2 



V, 



N 



and 



*(r)= [<£ ps (r) &(r) 4> 2 {v) ... ^(r) 



(6) 



(7) 



where A T denotes the transpose of A and <fr(i*o) = 
x P T (ro)V. Finally, we define the 12-component opera- 
tor 



V 



dx 



d_ 



dxdy dxdz 



dydx 



,)ij 2 



_ d 

Jz dzdx 



o 2 

dzdy 



t)z- 



(8) 



where the first three components are the gradient and 
the next nine components are the Hessian operator. 

The nine position, frequency, and axis constraints de- 
fined by Eqs.[3]and[5]can be assembled into one equation: 



C 1 (7'®* T (ro))V = C 2 , 



(9) 



where Ci is a j x 12 matrix and C 2 is a j X 1 column 
vector, where j — 9 for this example. 

The position constraints in Eq. [3] can be reconstructed 
by using Ci to pick out the three gradient components of 
V and C 2 to set them to zero. The constraints in Eq. [5] 
can be treated in a similar manner. Thus, to encode the 
nine desired constraints, we use 



100 000000000 
010 000000000 
001 000000000 



Ci = 



000 100000000 

000 000010000 

000 000000001 

000 010000000 

000 001000000 

000 000001000 



(10) 



c 2 = 







(m/q)u* 
{m/q)u£ . (11) 
(m/q)ul 







Additional white space has been inserted in both equa- 
tions to aid the reader by separating the position, fre- 
quency, and principal axis constraints into groups in the 
vertical direction, as well as separating the gradient and 
Hessian components of Eq.[l0]in the horizontal direction. 

Once Ci and C 2 arc determined, Eq. ([9| can be solved 
for V by inverting Ci (V^^ T (ro)) , thus determining the 
control voltages that create the desired trapping poten- 
tial. This inversion may not be strictly possible, as is 
the case when the number of constraints does not equal 
the number of control potentials, leading to an over- or 
under-determined problem. Also, we are interested only 
in solutions where the magnitudes of all control volt- 
ages are smaller than a maximal voltage V,nax (for our 
apparatus, V niax — 10 V). To achieve this, we use a 
constrained least-squares optimization algorithm, as de- 
scribed in Ref. 1271. to calculate 



min \CAV® ^ T (r ))V- C 2 | 

|Vi|<y max 1 v ' 1 



(12) 



In cases where Eq. ^ is over-constrained, this method 
yields a "best-fit" V. When Eq. Q is under-constrained, 
as is usually the case for large trap arrays with many 
electrodes, it returns a null space in addition to V, which 
can be added to V to find multiple independent solutions. 

Nine constraints were used above, but many are unnec- 
essary. For QIP in a linear trap array, constraining the 
axial mode frequency and orientation is often sufficient. 
Parameters for the other two modes are less important 
and often achieve reasonable values without being con- 
strained, in which case they can be omitted from the 
constraint matrices. 

In addition to explicitly defined user constraints, there 
are implicit physical and geometric constraints that must 
be considered. As an example, take the three secular fre- 
quencies of the ion, u> x , Lu y , and uj z . These frequencies 
result from a hybrid potential that includes both pseu- 
dopotential and control potentials. The contributions 
from both potentials can be separated mathematically 
into components, u>tt,i and taj respectively, which add in 
quadrature to give the overall frequency: ojf = u>i 2 +uif {i . 
(An imaginary frequency component would imply anti- 
trapping, while a real component yields trapping.) The 
control electrodes produce a quasi-static electric field, 
which Laplace's equation requires to be divergenceless. 
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This places a physical constraint on the frequencies com- 

3 

ponents due to the control potential, namely Hif = 0. 

t=i 

Thus, Laplace's equation permits only certain combina- 
tions of the secular frequencies. For a linear Paul trap, 
where w r f, z = 0, the secular frequencies must obey 



(13) 



where u> T f is the pseudopotential radial trapping fre- 
quency. 

The trap geometry can place constraints on the trap- 
ping potentials, as well. For example, in traps where the 
geometry contains some symmetry, the potentials must 
preserve that symmetry. Care must be exercised to en- 
sure that user-defined constraints do not contradict phys- 
ical or geometry constraints, as this will invalidate the 
solution. 

Though we invoke only position, frequency, and ori- 
entation constraints here, other varieties of user-defined 
constraints can be easily included with this framework, 
and a more complete discussion of these constraints is 
presented in [55]. The constraints used to construct the 
waveforms in the X-junction array were as follows: 

1. The position of the potential minimum was con- 
strained in three directions to be at ro. 

2. One of the principal axes was constrained to lie 
along the trap axis (which involves two constraints 
on axes orientation). 

3. The ion axial frequency was constrained (usually to 
3.6 MHz). 

4. The voltages were constrained to be between ±10 V 
(to conform to the limits of the voltage supplies 
used in the experiment). 

This relatively sparse set of constraints tended to give 
good solutions at most locations considered. Item 4 is an 
inequality constraint that is easily implemented by use 
of the constrained least-squares method. 

When solving waveforms that transport across multi- 
ple zones, V can become discontinuous from step to step, 
especially when transitioning between sets of control elec- 
trodes. These jumps occur when an under-constrained 
problem (with null space rank > 0) has multiple linearly- 
independent solutions and the algorithm returns a differ- 
ent solution from one step to the next: during transport 
there will be some position at which it is suddenly easier 
to produce the desired potential using a new combination 
of electrodes. In principle, such jumps should not have 
adverse effects on the potential at the ion, as the poten- 
tials on both side of the jump fulfill the same constraints 
and should transition smoothly. However, since the po- 
tentials on the electrodes are filtered, we would expect 
the potential at the ion to experience a transitory jump 
during the transition. 



These solution jumps can be handled by various means. 
We used the constrained least-squares method to seed 
each new solution with the solution of the previous step 
while introducing a cost for deviating from the previous 
solution by replacing Eq. [12] with 



* T (r ))V-C 2 



min CMT^ 

|V,|<V max , \Vi-V iMst \<a 1 

for a positive constant a. This removes the need for it- 
eratively choosing weights to keep the voltages within 
bounds, as suggested in Ref. [25] ■ This forced the jump 
transition to be extended over multiple steps, rather than 
allowing a discontinuous jump. Another approach is to 
average the two V's on each side of the discontinuity, 
taking advantage of the linearity of the equations, to 
produce an intermediate solution that still satisfies the 
constraints |28) . Performing several steps of such aver- 
aging will smooth the jump. Alternately, trial and error 
can often be used to determine a set of constraints that 
does not produce a jump, but this can require significant 
effort. 



1. Transport timing 

If the spatial interval between steps in the waveform 
is small enough, the potential, once applied to the elec- 
trodes, will move smoothly from step to step [26 . The 
velocity of the potential well (and, thus, the ion) is con- 
trolled by the rate at which the waveform steps are 
updated on the electrodes. In our case, the control 
potentials were supplied by digital-to-analog convert- 
ers (DACs) that had a constant update rate i?DAC = 
480 kHz and the number of update steps was adjusted to 
change the velocity. 

Different velocity profiles have been considered for 
minimizing excitation while transporting [29, 30, . In this 
report, the ions were usually transported by use of a con- 
stant velocity with equally spaced waveform steps. This 
could potentially lead to the ion being 'kicked' as the ve- 
locity jumps at the beginning and end of the transport, 
resulting in motional excitation. However, these velocity 
jumps were smoothed by low-pass filters placed on the 



control potentials (see Sec. VB). A smoother 'sinusoidal' 



velocity profile was also tested but was abandoned af- 
ter observing no discernible difference in the amount of 
motional excitation by use of the different profiles. This 
suggests that both transport protocols were well within 
the adiabatic regime at the speeds used. 

Low-pass filtering (160 kHz corner in our case) can 
also potentially distort the waveforms when transporting 
quickly, placing an upper limit on the ion speed. How- 
ever, the practical speed limit was set by the combination 
of the maximum update rate of the digital-to-analog con- 
verters and the number of update points required to ac- 
curately produce a continuous harmonic potential in the 
region of the pseudopotential barrier. This limit was ex- 
perimentally determined for each waveform by adjusting 
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FIG. 6: (a) The waveform (as a function of position rather 
than time) used when transporting an ion from the exper- 
iment zone, located at z = —880 /im, to the center of the 
junction at z — (im. We plot voltage versus the z position 
of the minimum of the trapping potential during the trans- 
port. The locations of the electrodes near the junction are 
depicted, along with their electrode number, by the rectan- 
gles in the bottom of the figure. The region from —100 to 
/im is inside of the junction. The voltage traces are num- 
bered to show which electrode they correspond to. Electrodes 
8 and 9 remain near V and are omitted for clarity. In ad- 
dition, the potentials applied to the control electrodes on the 
bottom wafer are not displayed, as they are nearly identical to 
those applied on the top wafer, (b) A schematic of the trap, 
showing the range over which this waveform transported. 



the number of update points until minimum motional 
excitation was achieved. If faster DACs are available 
and distortion of the waveforms due to low-pass filter- 
ing is of concern, the waveform can be pre-compensated 
to account for these distortions and produce the desired 
waveform at the ion. 

The waveforms used to transport from £ to C are dis- 
played in Fig. [6] as a function of the position of the min- 
imum of the trapping potential (the ion's location). The 
potentials applied to the lower trap-wafer control elec- 
trodes (on opposite sides of the main channel) were nearly 
identical and are omitted for clarity. These waveforms 
could be run left to right to transport an ion 880 fim from 
£ to C, or they could be run in reverse. The waveforms 
that transported ions into the other two branches of the 
junction (to J- and V) were similar to this waveform due 
to the approximate symmetry of the trap. 

In a typical transport, the potential minimum was 



moved at a constant velocity, and there was a direct lin- 
ear relationship between the location of the minimum 
(horizontal axis of Fig. ^&)) and the time elapsed since 
the beginning of the transport. The typical transport du- 
ration for the waveforms in Fig. |6][a) was approximately 
165 /is, with 50 /is to cross the pseudopotential barrier. 

Some control potentials reached the ±10 V limit placed 
by use of the constrained least-squares method while 
traversing the pseudopotential barrier near the junction. 
Other control potentials had sharp and abrupt changes, 
which resulted from the constraint in Eq. [14] that pre- 
vented 'solution jumping' by defining how much a given 
waveform step can deviate from the previous step. In- 
stead of jumping, the voltages ramped linearly over sev- 
eral steps. Although these individual potentials were not 
smooth in time, they were continuous, which was suffi- 
cient to ensure that the overall potential experienced by 
the ion evolved smoothly. 

The axial frequency was chosen to be 3.6 MHz and 
was held constant during much of the transport starting 
at £ and moving toward C (Fig. ]7). The frequency was 
(adiabatically) linearly ramped to 4.2 MHz as the ion ap- 
proached the apex of the pseudopotential barrier, making 
the ion less susceptible to rf-noise heating of the secular 
motion (see Sec. VA|. The value 4.2 MHz was the max- 
imum axial frequency attainable at the apex due to the 
strong anti-confinement of the pseudopotential at that 
location and the ±10 V limit of the DACs providing the 
control potentials. The axial frequency then continued to 
increase as the ion descended the barrier, reaching a final 
value of 5.7 MHz at C. At this location, all control poten- 
tials were V and the pseudopotential provided all the 
trapping, resulting in near-degenerate 5.7 MHz confine- 
ment along the x and z directions, while the y secular 
frequency was 11.3 MHz. When transporting multiple 
ions in the same potential well, it would be preferable 
to break the frequency degeneracy at C to ensure well- 
defined axes for the ions. In practice, the motional exci- 
tation rates when moving pairs of ions were still relatively 
low despite the near degeneracy at C (see Table |H| . 



IV. TRANSPORT EXPERIMENTS 

The transport experiments were performed with 
9 Be + ions inside a vacuum system with a pressure 
of p < 5 x 10~ n Torr = 7 x 10~ 9 Pa. A 1.3 x 
10~ 3 T magnetic field was applied to split the Zee- 
man states, and the ions were optically pumped to 
the 2 S 1/2 \F = 2,m F = -2) state (henceforth |2,-2)). 
Manipulation of the 9 Be + ion motional and internal 
states used the techniques of Refs. (6] [31]. Two-photon 
stimulated-Raman transitions enabled coherent transi- 
tions between the qubit states |2, —2) and |1, —1) at fre- 
quency loq rs 27T x 1.28 GHz. In addition, by tuning the 
difference frequency of the Raman beams to ujo ± , 
it was possible to drive a blue(red)-sideband transition: 
|2, -2) \n) <-» |1, -1) | ra ± 1). Here |n) is a Fock state of a 
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FIG. 7: Predicted secular frequencies as a function of position 
corresponding to the waveform in Fig. [6] The axial frequency 
along z was constrained to be 3.6 MHz during the major- 
ity of the transport, while the radial frequencies were uncon- 
strained. As the ion ascended the pseudopotential barrier, 
the axial frequency linearly ramped up to 4.2 MHz. Beyond 
the apex of the barrier, a second linear ramp was applied to 
bring the frequency up to 5.7 MHz. As the ion approached 
the center of the junction, the x and z frequencies became 
nearly degenerate. 



a stationary ion (~ 0.5/hr). Having observed millions 
of successive round trips for all three types of transport, 
combining all losses implies a transport success proba- 
bility of greater than 0.999995[1S]. Since transport com- 
prised a small fraction of the total experiment duration, 
many of these losses likely occurred when the ion was not 
being transported. In one instance, more than 1,500,000 
consecutive E-C-E transports were performed with a sin- 
gle ion. 

Loss rates for transported ion pairs were again compa- 
rable to stationary pairs (~ 2 per hour). Absolute pair 
loss rates were higher than those for single ions, presum- 
ably due to multi-ion effects [51 [33J . 

V. EXCITATION OF THE SECULAR MOTION 

Excitation of the ion's motion during transport was 
attributed to two main mechanisms: one due to rf noise 
near O r f and the other due to excitation from the digital- 
to- analog converters (DACs). 



selected motional mode. Ground-state cooling was per- 
formed by use of a series of red-sideband pulses, followed 
by repeated optical pumping to |2, —2). State detection 
was performed using state-dependent resonance fluores- 
cence, where predominantly the |2, —2) state fluoresces. 

Each transport began by cooling an ion (or ion pair) 
in £ to the motional ground state. The ion was then 
transported into or through the junction and returned to 
E. Three transports patterns were used: E-C-E moved 
to C and back, while E-C-T-C-E and E-C-V-C-E moved 
to T and V, respectively, before returning to £ . The £- 
C-£ transport moved the ion 1.76 mm, while E-C-T-C-E 
and E-C-V-C-E moved the ion 3.52 mm and 2.84 mm, 
respectively. Once the ion returned to £, the motional 
excitation was determined by measuring the asymmetry 
in red- and blue-sideband Raman transitions [311 H2] ■ 

To determine the single-ion transport success rate for 
E-C-T-C-E transports, two sets of 10,000 consecutive 
transport experiments were performed [T], but with the 
imaging system focused on £ in the first set and on J- in 
the second. The first set verified that the ion successfully 
returned to £ every time. The second set verified that 
the ion always reached T at the intended time. Together, 
these sets of experiments imply the success rate for going 
to T and returning to £ exceeded 0.9999. The procedure 
was repeated for E-C-V-C-E , with the same result. The 
E-C-E transport can not be verified in the same manner 
because the bridges obscure the ion at C, but since the 
ion must transport through this location to reach J- and 
V, the reliability should be no worse. 

Ion lifetime, and thus transport success probabil- 
ity, was ultimately limited by ion loss resulting from 
background-gas collisions [B] . With this in mind, the ion 
loss rate during transport was not larger than that for 



A. rf-noise heating 

Consider a trapping rf electric field with an additional 
sideband term, 

E rf (r, t) = E (r) [cos £l rf i + £ N cos (Q rf ± u z )t] , (15) 

where C 1, and O r f ± u z is at one of the two axial 
motional sidebands of the ion. In 1 , it was shown that 
the two terms will beat at uj z to produce a force that 
can excite the ion's motion. If the second term is not 
coherent, but instead is broad-spectrum noise, this will 
lead to excitation of the axial motion at a rate of 



<SV N (^rf 







dz 

00 z) 



E 2 (z) 



+ 



sw(firf - w z ) 



(16) 



where 5y N (f2 r f ±<j z ) is the voltage-noise spectral density 
at either the upper or lower rf sideband, and V T { is the 
amplitude of the trapping rf potential being applied to 
the rf electrodes. Eq (z) is the axial rf electric field ampli- 
tude at the location of the ion. This heating mechanism 
is proportional to the slope of the pseudopotential and is 
significant only in places with a large slope, such as the 
pseudopotential barriers near the junction (but not in, 
for example, £). 

This heating mechanism was verified in Ref. [1] by 
measuring the heating rate at various locations along 
the pseudopotential barrier between £ and C, while 
spectrally-dense white noise (centered on the lower side- 
band, fi r f — lo z ) was injected onto the trap rf drive. Fig- 
ure [8] plots the ratio of measured heating rate to esti- 
mated injected Sy N and theoretical values of this ratio 
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FIG. 8: The ratio of heating rate n z to voltage noise spectral 
density SV N (f2rf — u z ) for various locations along the trap axis 
(C is located at /im) . This figure is reproduced from Ref . [T| . 
The theoretical prediction used a pseudopotential modeled 
from electrode geometry and is shown both with and without 
a scaling parameter (= 1.4). The simulated pseudopotential is 
overlaid in the background, in units of eV. Since heating was 
gradient dependent, we saw very little heating at the peak of 
the pseudopotential barrier, even though this was the point of 
maximum (axial) rf electric field and therefore maximum ax- 
ial rf micromotion. Nearly identical pseudopotential barriers 
were present on the other three legs of the junction. 



according to Eq. ( 16 1 based on simulations of trap poten- 



tials, for the ion held at several positions between £ and 
C. A plot with the theoretical values multiplied by a scal- 
ing factor {— 1.4) is also included. The deviation of the 
scaling factor from 1 is not unreasonable due to the dif- 
ficulty of accurately measuring a variety of experimental 
parameters. 

The motional excitation for full junction transports 
to C was observed to decrease as the ion speed was in- 
creased, which minimized the exposure to the rf noise 
while on a pseudopotential slope. This continued up to 
a maximum speed limit, due to the slow DACs, above 
which the other excitation mechanism (below) began to 
dominate. At the optimum speed, the ion spent only 
approximately 50 /is on each barrier (above 10% of the 
barrier height). 

Another approach to mitigate rf noise is to suppress 
the sideband noise with better filtering of the applied 
rf trapping potential. In Ref. pQ, the large rf potential 
(Vrf « 200 Vpcak at fi rf m 2ir x 83 MHz) was provided 
by a series of tank resonators, which suppressed noise at 
the motional sidebands (±3.6 MHz). The primary res- 
onator was a quarter-wave step-up resonator |34) with a 
loaded Q = 42 and corresponding bandwidth (FWHM) 
of 2 MHz. This resonator extended into the vacuum, with 
the trap attached to the voltage anti-node. A second half- 
wave resonator with Q = 145 was attached, in series, to 
the input of the primary resonator, with a 3 dB atten- 
uator in between to decouple the two resonators. The 
resonant frequencies of the two resonators were tuned to 



be equal. This network resulted in an estimated ambi- 
ent SV N (fi r f ± to z ) of -177 dBc at the ion. In the work 
reported here, the second resonator was replaced with 
a pair of half- wave tank resonating cavities. This filter 
pair provided more than 38 dB suppression at frequencies 
f2 r f ± 27r x 3.6 MHz (when not coupled to the primary 
resonator), an additional suppression of approximately 
10 dB over the half- wave filter used in Ref. pQ. <SV N at 
the ion was not re-measured with this new filter pair, but 
observed reductions in excitation during transport were 
consistent with a 10 dB drop in rf noise. 



B. DAC update noise 

Another primary source of motional excitation was at- 
tributed to the 16-bit, ±10 V DACs that supplied the 
waveform potentials to the electrodes. The DAC voltages 
were updated at a constant rate Rbac (< 500 kHz), re- 
sulting in Fourier components that could excite the ion's 
motion if 27r x Rbac = ^z/J, for any integer J. 

This effect was observed by first preparing the ion 
in the motional ground state at £ and then transport- 
ing toward C. Instead of proceeding all the way to C, 
the transport was stopped (at z — —300 /im) before 
the axial frequency began to ramp up. Thus, the lo- 
cal potential- well frequency remained approximately con- 
stant at uj z = 2ir x 3.6 MHz. The ion was then returned 
to £. A red-sideband Raman 7r-pulsc for n — to n = 1 
excitation was applied to determine if the ion remained 
in the ground state [31] [32] . If the ion was excited out of 
the ground state during transport, the Raman pulse had 
a certain probability to transfer the ion into the |1, — 1) 
state, which did not fluoresce during detection. If the 
ion remained in the ground state, this side-band pulse 
had no effect and the ion remained in the bright |2, — 2) 
state. Thus, fluorescence detection after the side-band 
pulse could distinguish an excited ion from a non-excited 
ion. 

This experiment was performed for various values of 
Rbac, an d the results are shown in Fig. [9j It was diffi- 
cult to extract the ion's exact motional state after the 
transport, but the correlations between the ion's mo- 
tional excitation and u> z corresponding to a harmonic of 
Rdac were evident. The energy gain exhibited a reso- 
nance at several values for Rbac = oj z /(27rJ) with J = 8 
to 14. When the number of update steps was increased, 
while the update rate was held constant (which resulted 
in an increased transport duration), the bandwidth of 
these resonances decreased, as expected from a coherent 
excitation. 

Use of an update rate that was incommensurate with 
the motional frequency {Rbac 7^ w z/(27rJ)) minimized 
this energy gain. However, increasing the transport speed 
(using the same update rate) required a reduction in 
the number of waveform steps, which caused the reso- 
nances to broaden. Minimizing the rf-noise heating re- 
quired fast transport, so at the speed that gave the low- 
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FIG. 9: Plot showing the the number of fluorescence photons 
detected in a duration of 200 /is following round-trip transport 
and a subsequent red sideband pulse, for various DAC update 
rates Rdac- Before transport, the ion was prepared in the 
motional ground state and then transported through a specific 
waveform, where ui z — 2n x 3.6 MHz was maintained during 
the entire transport. If the ion remained in the ground state 
after transport, the ion fluorescence was at its maximum value 
(approximately 8 photon counts detected), but when the ion 
became motionally excited, the fluorescence dropped. As can 
be seen, the motion was excited at specific update frequencies 
that correspond to -Rdac = w z / (2tvJ) for J = 8 to 14 (marked 
by the vertical red lines). 



est rf-noise excitation rates, the DAC heating resonances 
were so broad that they overlapped, and there was no 
achievable -Rdac that would not result in energy gain. 
The DAC heating effect was further compounded by the 
fact that the axial frequency was not constant during a 
full junction transport, making it impossible to achieve 
-Rdac w z/(27rJ) for any constant -Rdac- The update 
frequency -Rdac = 480 kHz appeared to be most favor- 
able and was used for the results here and in Ref. PQ. 

Faster DACs capable of -Rdac > w z /27r should signif- 
icantly suppress this motional excitation. Alternatively, 
aggressive filtering of the DAC output can combat this 
problem. The results in Ref. [1] used the RC filter net- 
work shown in Fig. |10[ a), which provided suppression 
by two orders-of-magnitude over the range of lo z /2tt val- 
ues used during transport (3.6 to 5.7 MHz), but was not 
sufficient to completely suppress the DAC heating. In- 
creasing the RC time constant would increase the filtering 
but would also slow down the rate at which the ion can 
be transported. Instead, these simple RC filters were re- 
placed with the approximately third-order Butterworth 
filter 



shown in Fig. 10 'b). (The output impedance 



of the DAC was < 0.1 S2 and contributed minimally to 
the filter response.) 

A Butterworth filter has a frequency response given by 



1 



1 



\B n (iu/u )\ y/l + (uj/u} ) 2n ' 



(17) 




where B n (s) is the n th -order Butterworth polynomial and 
ll>o is the corner frequency. If n = 1, the frequency re- 



FIG. 10: The control voltages were provided by 40 indepen- 
dent DACs (only one shown here). The DAC output was 
filtered prior to being applied to the trap electrodes, through 
a two stage filter seen in (a). The control potentials were ref- 
erenced to the grounded vacuum system, which served as the 
rf ground as well, (b) After DAC-update noise was observed 
to excite the secular motion of the ions, the external filters 
were replaced by the approximate 3 rd -order 'Butterworth' fil- 
ter shown here. 



sponse reduces to a RC frequency response. In the ex- 
periments here, such higher-order filters provide stronger 
noise suppression at uj z while still allowing fast transport. 
A comparison of the theoretical response functions for the 
RC filters used in Ref. [1] and the Butterworth filters used 
here is shown in Fig. |11| The internal vacuum RC compo- 
nents already on the filter board were taken into account 
when planning the Butterworth filter, but the external 
filter components were designed to dominate the filter's 
response in the frequency range of concern. Thus, the fil- 
ter was approximately third-order, despite the presence 
of four components (including the filter board capacitor) 
with frequency-dependent impedances. The new filters 
increased the noise filtering by 22 dB at 3.6 MHz and 
26 dB at 5.7 MHz. Furthermore, the electric-field noise 
at the ion due to Johnson noise in the resistive elements of 
the new filters was less than that for the previous filters 
for all frequencies of interest. For the transport dura- 
tions used, these filters did not appreciably distort the 
waveform. 



C. Anomalous noise heating 

The 'anomalous heating', which is thought to arise 
from noisy electric potentials on the surface of the 
trap [32], was measured to be 40 quanta/s for ui z /2ir — 
3.6 MHz at £ and was not a significant source of ex- 
citation during the transport experiments. For exam- 
ple, we estimate that it should have contributed only 
0.007 quanta for £-C-£ transport. To compare the mea- 
surements of various ion traps, it has been common to 
express the heating in terms of the electric field noise 
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FIG. 11: Theoretical transfer function G{ui) versus fre- 
quency for the original RC filter (dashed) and the improved 
approximately-'Butterworth' filter (solid), where the new fil- 
ter had a faster roll-off at high frequency. Since the pass-band 
extended farther for the new filter, the transport speed could 
be increased while providing more filtering at the secular fre- 
quency (3.6 MHz to 5.7 MHz). Both traces include the RC 
components inside the vacuum. 

with the expression [52] : 



where S^(lj z ) is the spectral density of electric field fluc- 
tuations at the secular frequency. From the results here 
we found S e (2tt x 3.6 MHz) = 2.2 x 1(T 13 (V/m) 2 Hz _1 , 
where the distance of the ion to the nearest electrode sur- 
face was 160 //m. This result, when compared to other 
traps as in Refs. [TSJ l3l)H3"5] , was significantly below that 
of most other room-temperature ion traps. The cause of 
this relatively low heating rate is not known, but surface 
preparation could be a contributing factor. 

D. Other heating mechanisms 

In the experiments here, the transport was slow and 
the trapping potential changed slowly compared to the 
motional frequencies, so we did not expect non-adiabatic 
excitation of the motion. This was supported by observa- 
tions that the excitation did not decrease as the transport 
was slowed. For very slow transport, the heating actually 
increased because the ion spent more time crossing the 
rf barriers, resulting in increased rf-noise heating. Fur- 
thermore, no reduction in heating was observed when a 
gradual (sinusoidal) velocity profile was used instead of 
a constant velocity over the entire transport. 

The waveforms were produced assuming a specific 
value of V r { and corresponding pseudopotential. In the- 
ory, if the actual V r { does not match the assumed Vrf, 
the axial trapping potential will not be as intended at 
the barriers. In practice, there was an optimal value for 
the rf power which resulted in the lowest excitation rates 
and likely corresponded to the assumed V r [ . The rf power 
was prone to slow drifts over many minutes (likely due 
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S-C-S 


350 


0.053 ± 0.003 


310 


3.2 ± 1.8 


E-C-T-C-E 


910 


0.18 ±0.02 


630 


7.9 ± 1.5 


E-C-V-C-S 


950 


0.18 ±0.02 


870 


14.5 ±2.0 



TABLE I: The axial-motion excitation An for a single 
9 Be + ion for three different transports through the X- 
junction. The results of this work, as well as that of Ref. PQ, 
are given for comparison^] The transport duration includes 
20 /xs for the ion to remain stationary at the intermediate des- 
tination (30 fis for the data from Ref. pQ), before returning to 
£ . The energy gain per trip is stated in units of quanta in a 
3.6 MHz trapping well where An = 0.1 quantum corresponds 
to 1.6 neV. 

"The transport durations given in Ref. [Tj were reported in error. 
The correct values are 140 fis for transporting from £ to C, 300 fis 
to go from £ to V, and 420 fis to go from £ to T . This error did 
not affect any other results in Ref. pp. 

to temperature drifts in the resonators) which resulted 
in modest increases in motional excitation; it was nec- 
essary to occasionally adjust the rf power (every 10 to 
30 minutes) and hold it to within < 1 % to achieve the 
lowest motional-excitation rates. In practice, this was 
performed by ensuring that the radial secular frequen- 
cies at £ remained constant. 



E. Motional excitation rates 

The motional excitation for single-ion transports was 
measured by use of sideband asymmetry measure- 
ments [31] after a single pass through the junction, and 
the results are summarized in Table ??. These results 
were significantly better than those in Ref. PQ , which are 
listed for comparison. In Ref. [1], rf noise was estimated 
to contribute 0.1 to 0.5 quanta of excitation per pass over 
a pseudopotential barrier, which explained between 3 and 
30% of the excitation seen. The remainder of the excita- 
tion was attributed primarily to DAC update noise. The 
new trap rf filters and control electrode Butterworth fil- 
ters produced the observed reduction in excitation rates. 

The transport durations, which were optimized for 
minimal excitation, are also given in Table ??. The tabu- 
lated durations correspond to the full transport duration 
including returning to £ (rather than the half-transport 
reported in Ref. pp. The durations also include a 20 fis 
wait at the half-way point (C, J- , or V, depending on the 
transport) for the new results and a 30 /is wait for those 
from Ref. PQ. 

Moving pairs of ions in the same trapping well would 
be useful for both sympathetic cooling and efficient ion 
manipulation [7J. This type of transport was demon- 
strated by use of pairs of 9 Be + ions and the measured 
motional excitation is reported in Table [TTJ Excitation in 
both the center-of-mass (COM) and stretch modes was 
measured. Additional heating mechanisms for multiple 
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COM 


Stretch 


COM 


E-C-E 


0.39 ± 0.03 


0.13 ±0.02 


5.4 ± 1.2 


E-C-T-C-E 


0.67 ±0.05 


0.53 ±0.05 


16.6 ± 1.8 


E-C-V-C-E 


0.72 ±0.06 


0.14 ± 0.02 


53.0 ± 1.2 



TABLE II: The axial-motion excitation An for a pair of 
9 Be + ions transported in the same trapping well. Values for 
both axial modes of motion (COM and stretch) are reported. 
The energy gain per trip is stated in units of quanta where 
the COM frequency is 3.6 MHz and the stretch frequency is 
6.2 MHz. Results from Ref. [T] are also given, though only 
the COM mode excitation was investigated. 

ions [SI may explain the higher energy gain observed 
for the pair. For E-C-V-C-E transport, the two-ion crystal 
must rotate from the z axis to the x axis and back. For 
the waveforms used, the potential was nearly the same in 
the x and z directions at C. Therefore the axes were not 
well defined throughout the transport, which can lead to 
an uncontrolled rotation of axes. It is possible that the 
discrepancy in the excitation between E-C-T-C-E and E- 
C-V-C-E for two ions may have resulted from this uncon- 
trolled rotation at C. 

We expect (and observed) less excitation of the stretch 
mode relative to the COM mode, for two reasons. First, 
the stretch mode frequency was higher than that of the 
COM mode (wstr = V3wcom)- Thus, the filters on the 
rf and control potentials were more effective at suppress- 
ing noise that could excite the stretch mode. Second, a 
stretch mode can be excited only by a differential force on 
the two ions, while the COM mode is excited by a force 
common to both ions. Given the proximity of the ions to 
each other (a few micrometers) compared to the distance 
of the ions to the trap electrodes, the relative amplitude 
of differential forces acting on the ions are expected to 
be less than common forces. 



VI. MITIGATION OF MAGNETIC FIELD 
FLUCTUATIONS 

So far, we have discussed the suppression of undesired 
excitation of motional degrees of freedom. We now dis- 
cuss how magnetic-field fluctuations affecting internal- 
state (qubit) coherence are suppressed in the X-junction 
trap array. 

Decoherence of superpositions of the \2, — 2) and 
1 1,-1) qubit basis states occurs both during transport 
and while the qubit is stationary. Previous experiments 
demonstrated that junction transport contributed negli- 
gibly to decoherence [T] . Magnetic field fluctuations form 
the dominant contribution to qubit dephasing, yield- 
ing typical values (in this trap and others) of less than 
100 /is [39] . Use of a magnetic-field-insensitive qubit con- 
figuration enables extension of the coherence time to ap- 
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FIG. 12: Mu-metal magnetic shield. The main dome en- 
closed both the trap and the magnetic-field coils. Cylindrical 
tubulation extended along a glass vacuum envelope, which 
corresponds to the z direction at the trap. Reentrant flanges 
minimized field leakage around the imaging and optical access 
points. 

proximately 10 s and can be used with some gate opera- 
tions such as the M0lmer-S0renson gate [40] , but excludes 
implementation of a z a z gates [3"9l |4"T1 14"2"] . 

To suppress the effects of magnetic field fluctua- 
tions, we enclosed the trap and field coils inside a 
high-magnetic-susceptibility mu-metal shield and imple- 
mented an active magnetic field stabilization system. 
The shield (Fig. |l~2"| ) was designed for compatibility with 
the existing trap vacuum envelope and optical systems, 
and for ease of installation without the need to lift the 
trap apparatus from the supporting table. A cylindri- 
cal body and approximately hemispherical dome were 
selected based on general guidelines for magnetic shield- 
ing and manufacturing constraints. The main body 
and baseplate of the shield were constructed of 3.2 mm 
thick, single-layer mu-metal in order to provide maximum 
shielding of low- frequency magnetic field fluctuations and 
to suppress magnetic saturation of the mu-metal. This 
latter constraint arose because part of the field coils defin- 
ing the quantization axis of the qubits were located ap- 
proximately 1 cm from the walls, where a calculated field 
of 4 x 10~ 3 T was expected for typical operating condi- 
tions. 

Openings in the shield for optical components or laser 
beams were outfitted with a reentrant flange fastened to 
the main body of the shield. The flanges extended both 
outwards and inwards in the shield in order to maximize 
flux-line redirection. A target length-to-diameter ratio 
of 5:1 guided design but was typically not achieved due 
to geometric constraints arising from the exterior dimen- 
sions of the vacuum envelope and the desire to position 
optical elements as close as possible to the shield for max- 
imum beamline stability. These flanges were designed to 
be modular, allowing for redesign and replacement if in- 
creased shielding became necessary. Magnetic continuity 
was achieved for all mating flanges by use of internally 
threaded fasteners, producing a snug contact fit. 

The structure was designed to provide a minimum 
22 dB shielding of low-frequency fields. This was con- 
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firmed by use of a pickup coil and detecting 60 Hz fluc- 
tuations. The minimum suppression measured for the 
shield alone was > 20 dB parallel to the z axis of the 
trap, and the direction of the largest access opening in 
the shield. Shielding in excess of 60 dB at 60 Hz was 
measured in transverse directions. 

The 1.3 x 10~ 3 T quantization field is oriented 45° 
with respect to the vacuum envelope axis. Initial mea- 
surements identified current instability due to power- 
supply ripple as dominating measured decoherence once 
the shield was installed. We implemented a cus- 
tom current-regulation system based on a proportional- 
integral-diffcrcntial feedback circuit, a series current- 
sense resistor, and a low-current field-effect transistor. 
To minimize the effect of thermal drifts in the electronic 
circuit on magnetic field stability, we selected special low- 
TC (thermal coefficient) components and temperature- 
stabilized the enclosure. The most critical components 
were the gain and sense resistors; these were selected to 
be low-TC metal foil resistors with less than 2 ppm/K 
and less than 3 ppm/K stability, respectively. A four- 
terminal current-sense resistor was selected with high- 
power-handling construction (0.1 for the main coil and 
0.25 fl for the transverse shim coils). Similar care was 
taken to select low-TC difference amplifiers for the in- 
put stage and a low-TC voltage reference. All sense and 
feedback components were thermally sunk to an Al en- 
closure that was thermally stabilized by use of Peltier 
coolers and a commercial temperature controller with 
milliKelvin stability. Stabilization reduced current ripple 
from ~ 1 mA to ~ 30 /iA on the main field-coil current 
of 1.2 A. Net magnetic field fluctuations due to current 
ripple at the location of the ions were ~ 26 nT. 

Measurements of the dephasing time including both 
the magnetic shield and the stabilization circuitry 
demonstrated extension of the qubit coherence to 1.41 ± 
0.09 ms, more than 15 x longer than that without shield- 
ing and current stabilization, and sufficient for multiple 
transports before the qubit dephased. A spin-echo pulse 
doubles the coherence time to 2.99 ± 0.04 ms, indicating 
that slow shot-to-shot field fluctuations are small, and 
that decoherence is dominated by fluctuations on a mil- 
lisecond time scale. 



VII. MODE ENERGY EXCHANGE 

The secular modes of the ions were constrained to 
change throughout the transports, both in frequency 
and orientation. For most parts of the transport, the 
splittings between the mode frequencies were sufficiently 
large and the transport speed was sufficiently slow that 
modes changed adiabatically and energy did not trans- 
fer between modes. However, at C, the two princi- 
ple axes that lie in the (x, z) plane were designed to 
have nearly degenerate secular frequencies (uj' x ~ uj' z ~ 
2ir x 5.7 MHz), which could lead to mode-mixing. The 
third mode along y had a significantly higher frequency 



(jj' y = 2tt x 11.3 MHz, and would remain decoupled from 
the x and z modes. Since the radial modes were only 
Doppler laser-cooled before transport, x/z mode mix- 
ing would increase the excitation of the axial mode dur- 
ing transport. We employed two approaches that would 
minimize such axial excitation. First, the duration dur- 
ing which the ion was at C could be adjusted such to 
minimize the energy transfer between modes (by using 
a duration that corresponded to a full cycle of the mix- 
ing process). Alternately, a potential could be applied 
to various electrodes, which we will call the shim po- 
tential, to sufficiently break the degeneracy (in practice, 
\uj' x — ui' z \ > 27rx400 kHz could be achieved) and suppress 
the mixing. Both methods were effective and yielded 
similar transport excitation, though the second approach 
was used for the results in Tables ?? andlTTl 

However, in separate experiments, we explored a 
method for controlling energy transfer between the mo- 
tional modes of a single ion by using field shims near the 
junction to tune oj' x and w' z to near-degeneracy. Ideally, 
a demonstration of the method would work as follows. 
Prior to transport from £, the ion is cooled to the axial 
ground state \n z =0) along z and prepared in Doppler- 
cooled thermal states in the transverse modes. If the rel- 
ative orientation of the modes remains stationary as the 
ion approaches C, the modes should not exchange energy, 
even if they become degenerate. However, if the x and 
z mode directions diabatically (fast compared to 1/Acj) 
rotate 45° to new directions given by x' = (x+z)/y/2 and 
z' = (x — z)/^/2, we would expect the initial x oscillation 
to project onto the new mode basis with half of the energy 
going into each of the new modes. If Auj' = oj' x — uj' z ^ 0, 
the two oscillations would then begin acquiring a relative 
phase cj) = Auj't, where t is the period spent at C. By 
then quickly transporting away from C towards £ such 
that the mode axes rotate diabatically by —45° back to 
their original orientation, the oscillations would project 
back onto the original oscillator basis. If the wait period 
is such that <f> = ir X M (where M is an integer), the 
motion originally in the x mode would project back into 
the same mode. If, however, <f> — ^ x (2M — 1), then 
the ir-motion would project into the z mode; that is, the 
energy would exchange between x and z modes. 

We demonstrated the basic features of this exchange 
as follows. To tune the ui' x and uj' z close to degeneracy, an 
external shim potential with an adjustable amplitude was 
applied. The shim potential consisted of various contri- 
butions from 17 control electrodes near C, each multiplied 
by the overall scaling factor A. These contributions were 
selected so that the net shim potential would primarily 
alter the frequency splitting without significantly affect- 
ing other trapping parameters (such as the position of 
the trapping minimum and the y mode frequency). 

The black trace in Fig. [l3|^a) shows n z for the ax- 
ial mode after £-C-£ transport, as the wait period at C 
was varied. We derive n z from sideband measurements, 
as described above, and assume a thermal distribution. 
Although this assumption may not be strictly valid, it 
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FIG. 13: (a) Average motional excitation (n z ) in the axial 
mode after an £-C-£ transport, versus the duration at C (wait 
period). The black trace indicates exchange of energy be- 
tween the z mode, prepared near the ground state (at £), and 
the x mode, prepared in a thermal state via Doppler cooling. 
The smaller blue trace represents the identical preparation, 
except the transport was performed twice. During the first 
transport, the wait period was set to maximize the energy 
transfer from the radial mode to the axial mode, followed by 
returning the ion to £. After re-cooling the axial mode to 
the ground state, the round-trip transfer was repeated. The 
contrast was decreased (blue trace), indicating that less trans- 
verse mode energy was available for transfer to the z mode 
and therefore indicating cooling of the x mode, (b) and (c) 
The exchange contrast and exchange frequency (respectively), 
plotted versus the shim-potential scaling factor A. 



should give a reasonable approximation for n z < 1. Os- 
cillations between an excited and near-ground state en- 
ergy are visible. The projection process began prior to 
the ion reaching C and transport was too slow for the 
projection to be perfectly diabatic. Thus, the phase of 
the exchange oscillation in Fig. 
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is not well deter- 
mined and was observed to depend on both the exchange 
frequency and the details of the approach to C (includ- 
ing speed and trajectory). In practice, it was difficult 
to maintain a constant phase for more than a few min- 
utes, as drifts in the potential, likely caused by transient 
charge buildup and dissipation on the electrodes and also 
pscudopotential amplitude changes, caused the exchange 
frequency to drift over that time scale. 



Figure 13 'b) displays the oscillation contrast An z = 
max(n z ) — min(ri z ) for various scaling factors, A, of 
the shim potential, while Fig. |13[ c) gives the frequency 
of those oscillations versus the shim scaling factor. In 
separate experiments, the two mode frequencies at C 
were measured as a function of A by driving excitations 
with an oscillatory potential applied to the control elec- 
trodes. The difference between the two mode frequencies, 
Aw' = \ui' x — uj' z \, was observed to match the oscillation 
frequency of the exchange process. Figure 13 ^c) suggests 
Aw' is high on the extreme ends of the A range, while 
Fig. 13 [b) shows a reduction in contrast in these regions 
of high Aw', likely due to the reduction in diabaticity 
when Aw' was large. This conclusion was supported by 
the observation that the contrast decreased as the ion 



transport speed was reduced. However, there was a max- 
imum speed, above which contrast no longer increased, 
because other sources of excitation began to obscure the 
oscillatory signal. 

From Fig. 13 b), we see that the fringe contrast was 
also minimized for shim scaling factors near A — —0.15, 
where Aw' was small. This reduction in contrast can 
be explained as coinciding with the condition where the 
initial mode orientation is identical to the rotated mode 
orientation and thus the modes do not mix when pro- 
jected, which is a condition not necessarily related to 
Aw'. (We note that the A value for minimum exchange 
frequency does not match that for minimum contrast.) 
Unfortunately, it was not possible to verify this, as we 
could not measure the mode orientation at C, due to lack 
of laser-beam access. 

In the case where the energy from the x mode was 
transferred into the z mode, the ion could be returned 
to £ for a second round of ground-state cooling of the z 
mode. In the experiment, the exchange process in C was 
repeated and the results are shown as the blue trace in 
Fig. [13^ a), where a noticeable decrease in the ion's axial 
excitation was observed compared to the first experiment 
without the second stage of cooling. The small relative 
phase shift for the two traces in Fig. [l3^a) was due to 
the slow drift of Aw' over the several minutes required 
to take the two traces. 

Ideally, all of the energy would be transferred from the 
x mode into the z mode, and the subsequent cooling of 
the z mode would leave both modes in the ground state, 
leading to no oscillation during the second trip into C. 
In practice, complete transfer was inhibited for two pri- 
mary reasons. First, the ions were not being transported 
fast enough to make a clean diabatic projection of the 
motion onto the switched axes. Second, for complete en- 
ergy transfer, the projection should be onto axes rotated 
by ±45°. Any other angles would have resulted in in- 
complete transfer of energy. Attempts were made to ad- 
just additional shims in hopes of realizing configurations 
closer to ±45°. However, as Aw' — > 0, it, again, becomes 
difficult to predict the mode orientation with our ideal- 
ized computer models and we could not experimentally 
determine the mode axes in C. 

Nevertheless, we observed a clear and easily repro- 
ducible reduction in maximum oscillation amplitude 
(Fig. [l3^a)) from max^) = 0.68 ± 0.08 to max(n 2 ) = 
0.40 ± 0.05 after the second round of cooling, indicating 
the radial mode energy was being reduced. The use of ad- 
ditional rounds of exchange followed by cooling reduced 
max(n z ) further, but after three or four exchange rounds, 
other sources of excitation offset the energy reduction. 

When optimized, this technique might be used to cool 
all modes of a single ion to the ground state, while having 
the ability only to ground-state cool a single mode, as for 
the laser beam configuration used here. A junction is not 
required; all that is needed is a trap that can diabatically 
change the relevant mode orientations by ±45°, which 
could be possible in many trap configurations. 
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VIII. CONCLUSION 

In conclusion, we have demonstrated that transport 
through a two-dimensional trap array incorporating a 
junction can be highly reliable and excite an ion's motion 
by less than one quantum. This is a significant improve- 
ment over prior work with junction arrays (HUH] and sug- 
gests the viability of trap arrays incorporating junctions 
for use in large-scale ion-based quantum information pro- 
cessing. In addition, we have implemented a mu-metal 
shield and current stabilization to reduce qubit decoher- 
ence. We have also examined a technique for transferring 
energy between motional modes. 

DAC-update noise can be mitigated with the use of 
more appropriate filters such as the Butterworth filters 
used here and/or faster DAC update rates. Noise on the 
trapping rf potential can result in motional excitation at 



pseudopotential barriers as described in Sec. V A The 
junction design criteria in Ref. j!9j included minimizing 
these barriers. However, the results here show that the 
slope of the pseudopotential barrier is more important 
than the barrier height for suppressing motional excita- 
tion, suggesting that suppression of barrier height may 
not be a necessary constraint in future designs. Also, 
as observed here, with proper rf filtering, significant bar- 
rier slopes can be tolerated without causing significant 
heating. 

The technique for determining the waveforms de- 
scribed in this report can be extended to incorporate 
multiple trapping wells by expanding Eq. ^ to include 
multiple minima. Transport procedures such as the ion 



exchange in Ref. [T7j are also amenable to these solving 
techniques. Separating and combining of trapping wells 
requires consideration of the potential's quartic term |43j . 
Therefore, V in Eq. {[sj) could be expanded to include 
fourth-order derivatives. 

With the use of multiple junctions, the techniques de- 
scribed here could help provide a path toward transfer of 
information in a large-scale ion-based quantum processor 
and enable an increased number of qubits in quantum al- 
gorithm experiments. To do this, waveforms must be 
expanded to incorporate many trapping wells. Also, it 
is likely that a sympathetic cooling ion species will need 
to be co-trapped with the qubit ions to allow removal 
of the motional excitation from electronic noise, multi- 
ple junction transports, and separating and re-combining 
wells [HI EH HI] . If sympathetic-cooling ions are present, 
it may be advantageous to transport both ion species 
through a junction in a single local trapping well. Since 
the pseudopotential (and micromotion) are mass depen- 
dent, the qubit and cooling ions will experience different 
potentials, which could lead to additional motional ex- 
citation. If such excitation is excessive, it should still 
be possible to separate the ions into individual wells by 
species and pass the different species through junctions 
separately, followed by recombination. 

We thank J.M. Amini, K.R. Brown, and J. Britton 
for contributions to the apparatus and R. Bowler and 
J. P. Gaebler for comments on the manuscript. This work 
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Quantum Information program. This is a publication of 
NIST and not subject to U.S. copyright. 



[1] R.B. Blakestad, C. Ospelkaus, A. P. VanDevender, 
J.M. Amini, J. Britton, D. Leibfried, and D.J. Wineland, 
Phys. Rev. Lett. 102, 153002 (2009). 

[2] R. Blatt and D.J. Wineland, Nature 453, 1008 (2008). 

[3] C. Monroe and M.D. Lukin, Physics World 21, 32 (2008). 

[4] H. Haffner, C.F. Roos, and R. Blatt, Phys. Rep. 469, 
155 (2008). 

[5] J.I. Cirac and P. Zoller, Nature 404, 579 (2000). 

[6] D.J. Wineland, C. Monroe, W.M. Itano, D. Leibfried, 

B.E. King, and D.M. Meekhof, J. Res. Natl. Inst. Stand. 

Technol. 103, 259 (1998). 
[7] D. Kielpinski, C. Monroe, and D.J. Wineland, Nature 

417, 709 (2002). 
[8] L.-M. Duan and C. Monroe, Rev. Mod. Phys. 82, 1209 

(2010). 

[9] G.-D. Lin, S.-L. Zhu, R. Islam, K. Kim, M.-S. Chang, 
S. Korenblit, C. Monroe, and L.-M. Duan, Euro. Phys. 
Lett. 86, 60004 (2009). 

[10] J.I. Cirac and P. Zoller, Phys. Rev. Lett. 74, 4091 (1995). 

[11] R.J. Hughes, D.F.V. James, E.H. Knill, R. Laflamme, 
and A.G. Petschek, Phys. Rev. Lett. 77, 3240 (1996). 

[12] M.A. Rowe, A. Ben-Kish, B. Demarco, D. Leibfried, 
V. Meyer, J. Beall, J. Britton, J. Hughes, W.M. Itano, 
B. Jelenkovic, C. Langer, T. Rosenband, and 
D.J. Wineland, Quant. Inf. Comput. 2, 257 (2002). 



[13] M.D. Barrett, J. Chiaverini, T. Schaetz, J. Britton, W.M. 

Itano, J.D. Jost, E. Knill, C. Langer, D. Leibfried, R. Oz- 

eri, and D.J. Wineland, Nature 429, 737 (2004). 
[14] G. Huber, T. Deuschle, W. Schnitzler, R. Reichle, 

K. Singer, and F. Schmidt-Kaler, New J. Phys. 10, 

013004 (2008). 

[15] J. P. Home, D. Hanneke, J.D. Jost, J.M. Amini, 
D. Leibfried, and D.J. Wineland, Science 325, 1227 

(2009) . 

[16] J.F. Eble, S. Ulm, P. Zahariev, F. Schmidt-Kaler, and 
K. Singer, J. Opt. Soc. Am. 27, A99 (2010). 

[17] F. Splatt, M. Harlander, M. Brownnutt, F. Zahringer, 
R. Blatt, and W. Hansel, New J. Phys. 11, 103008 (2009). 

[18] W.K. Hensinger, S. Olmschenk, D. Stick, D. Hucul, 
M. Yeo, M. Acton, L. Deslauriers, C. Monroe, and 
J. Rabchuk, Appl. Phys. Lett. 88, 034101 (2006). 

[19] J.M. Amini, H. Uys, J.H. Wesenberg, S. Seidelin, J. Brit- 
ton, J.J. Bollinger, D. Leibfried, C. Ospelkaus, A. P. Van- 
Devender, and D.J. Wineland, New J. Phys. 12, 033031 

(2010) . 

[20] D.L. Moehring, C. Highstrete, M.G. Blain, K. Fortier, 
R. Haiti, D. Stick, and C. Tigges, arXiv: 1105. 1834 (2011). 

[21] J.D. Jost, Ph.D. thesis, University of Colorado (2010). 

[22] J. Chiaverini, R.B. Blakestad, J. Britton, J.D. Jost, 
C. Langer, D. Leibfried, R. Ozeri, and D.J. Wineland, 



16 



Quant. Inf. Comput. 5, 419 (2005). 
J.H. Wesenberg, Phys. Rev. A 79, 013416 (2009). 
M.N.O. Sadiku, Numerical Techniques in Electromagnet- 
ics (CRC Press, Boca Raton, 2009), 3rd ed. [38] 
K. Singer, U. Poschinger, M. Murphy, P. Ivanov, 
F. Ziesel, T. Calarco, and F. Schmidt-Kaler, Rev. Mod. 
Phys. 82, 2609 (2010). [39] 
R.B. Blakestad, Ph.D. thesis, University of Colorado 
(2010). 

P.E. Gill, W. Murray, and M.H. Wright, Practical Opti- 
mization (Academic Press, London, 1981). 

H. Uys, private communication. [40] 

R. Reichle, D. Leibfried, R.B. Blakestad, J. Britton, 

J.D. Jost, E. Knill, C. Langer, R. Ozeri, S. Seidelin, and [41] 

D.J. Wineland, Fortschr. Physik 54, 666 (2006). 

D. Hucul, M. Yeo, S. Olmschenk, C. Monroe, W.K. 

Hensinger, and J. Rabchuk, Quant. Inf. Comput. 8, 501 

(2008). [42] 

C. Monroe, D.M. Meekhof, B.E. King, S.R. Jefferts, 
W.M. Itano, D.J. Wineland, and P. Gould, Phys. Rev. 

Lett. 75, 4011 (1995). [43] 
Q.A. Turchette, D. Kielpinski, B.E. King, D. Leibfried, 

D. M. Meekhof, C.J. Myatt, M.A. Rowe, C.A. Sackett, [44] 
C.S. Wood, W.M. Itano, C. Monroe, and D.J. Wineland, 

Phys. Rev. A 61, 063418 (2000). 

H. Walther, Adv. At. Mol. Phys. 31, 137 (1993). [45] 
S.R. Jefferts, C. Monroe, E.W. Bell, and D.J. Wineland, 
Phys. Rev. A 51, 3112 (1995). 

U. Tietze and Ch. Schenk, Electronic Circuits: Handbook 
for Design and Applications (Springer, Berlin, 2008). 
L. Deslauriers, P.C. Haijan, P.J. Lee, K.A. Brickman, 
B.B. Blinov, M.J. Madsen, and C. Monroe, Phys. Rev. 
A 70, 043408 (2004). 

R.J. Epstein, S. Seidelin, D. Leibfried, J.H. Wesenberg, 
J.J. Bollinger, J.M. Amini, R.B. Blakestad, J. Britton, 



J. P. Home, W.M. Itano, J.D. Jost, E. Knill, C. Langer, 
R. Ozeri, N. Shiga, and D.J. Wineland, Phys. Rev. A 76, 
033411 (2007). 

N. Daniilidis, S. Narayanan, S.A. M0ller, R. Clark, 
T.E. Lee, P.J. Leek, A. Wallraff, St. Schulz, F. Schmidt- 
Kaler, and H. Haffner, New J. Physics 13, 013032 (2011). 

C. Langer, R. Ozeri, J.D. Jost, J. Chiaverini, B. De- 
Marco, A. Ben-Kish, R.B. Blakestad, J. Britton, 

D. B. Hume, W.M. Itano, D. Leibfried, R. Re- 
ichle, T. Rosenband, T. Schaetz, P.O. Schmidt, and 
D.J. Wineland, Phys. Rev. Lett. 95 (2005). 

K. M0lmer and A. S0rensen, Phys. Rev. Lett. 82, 1835 
(1999). 

D. Leibfried, B. DeMarco, V. Meyer, D. Lucas, M. Bar- 
rett, J. Britton, W.M. Itano, B. Jelenkovic, C. Langer, 
T. Rosenband, and D.J. Wineland, Nature 422, 412 
(2003). 

P.J. Lee, K.-A. Brickman, L. Deslauriers, P.C. Haijan, L.- 
M. Duan, and C. Monroe, J. Opt. B: Quantum Semiclass 
Opt. 7, S371 (2005). 

J. P. Home and A.M. Steane, Quant. Inf. Comput. 6, 289 
(2006). 

D. Hanneke, J. P. Home, J.D. Jost, J.M. Amini, 
D. Leibfried, and D.J. Wineland, Nature Phys. 6, 13 
(2009). 

For the 0.999995 success probability figure, we only ver- 
ified that each transport successfully returned the ion to 
£ , not whether the ion successfully moved the ion through 
the junction. Given the low rates of motional excitation, 
and the fact that we did verify that the ions move through 
the junction using 10,000 experiments, it is reasonable to 
assume the ion did travel through the junction if the ion 
successfully returned to £ . 



